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We investigate the effects of a finite chemical potential 
(i in QCD viewed as a disordered medium. In the quenched 
approximation, A4 = ip induces a complex electric Aharonov- 
Bohm effect that causes the diagonal contribution to the 
quark return probability to vanish at p = (half the pion 

mass). In two-color QCD, the weak-localization contribution 
to the quark return probability remains unaffected causing a 
mutation in the spectral statistics. In full QCD, the complex 
electric flux is screened and the light quarks are shown to 
diffuse asymmetrically with a substantial decrease in the con- 
ductivity along the 'spatial' directions. Mean-field arguments 
suggest that a d=l percolation transition may take place in 
the range 1.5po < p < 3po, where po is nuclear matter density. 

PACS numbers : 11.30.Rd, 12.38.Aw, 64.60. Cn 



1. QCD at finite chemical potential p is still not well 
understood despite the many efforts invested by a num- 
ber of groups in the past years. Lattice Monte-Carlo 
algorithms are difficult to implement at finite p owing 
to the complex character of the measure jl]]. Results 
from strong coupling arguments and quenched simu- 
lations [Qll are available but do not seem to be transpar- 
ent physically. A number of theoretical constraints can 
be implemented using symmetry and data B. However, 
they are only reliable for densities typically of the order 
of nuclear matter density pq. Results from constituent 
quark models at higher densities point at the possibility 
of a chiral transition at about 3 times nuclear matter den- 
sity j5| and the occurrence of a diquark superconducting 
phase at even higher densities ](|. 

In this letter we would like to address the effects of a fi- 
nite chemical potential on the chiral disorder of the QCD 
ground state. At p — we have recently [pj,B shown that 
light quarks in a finite Euclidean volume V are in a diffu- 
sive mode, with a diffusion constant D = 2F 2 /£ where 
F is the weak pion decay constant and S = | (qq) | the 
light quark condensate. The effects of matter cause the 
medium to change thereby affecting the diffusion prop- 
erties of the light quarks. In many ways our problem is 
similar to the problem of electrons in disordered metals 
under the influence of external sources || . 

2. The eigenvalue equation of the Dirac operator for 
fundamental quarks in a fixed gluon field A at finite 
chemical potential p is 

(if[A] + iptji) q k = Xk[A] q k • (1) 



for the right-eigenfunctions, and 

(if[A] - ^74) Qk = \* k [A]Q k . 



(2) 



for the left-eigenfunctions. The eigenvalues are com- 
plex and paired by chiral symmetry. The set (qk,Qk) 
is biorthogonal. Generalizing the construction Q for the 
case of a finite chemical potential, we may write the prob- 
ability p(t, p) for a light quark to start at x(0) in V and 
return back to the same position x(t) after a proper time 
duration t, as 

p(t,p) = ^e-^I^Kx^le^l+^l^xCO))! 2 ) . (3) 
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The averaging in ([5]) is over all gluon configurations us- 
ing the unquenched QCD measure with massive (sea) 
quarks. The normalization in @ is per state, where N 
is the total number of quark states in the four-volume V. 
Equation (0) can be resolved in terms of (0-0) 
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x / e <|t|(A, -K)W q . {x)Q * {x) Qk {x)ql {x] 



(4) 



where the exponent e _2m ''' is solely due to the valence 
quark mass. We note that (Q) is gauge-invariant and 
amenable to lattice Monte-Carlo simulation. It requires 
both the eigenvalues and eigenfunctions. 

For analytical considerations, it is best to rewrite (|^) 
in terms of the standard Euclidean propagators for the 
quark field, 

p(i, M ) = Z!lim [ ^^ e -i(^-x,)\t 



N y~xj (2tt) 2 

x^Tr (S(x,y; z u p)S\x,y; z 2 , p)) ) (5) 



with z\ 2 = m — i\\ 2, and 



S(x,y;z,fi) = (x\ 



if[A] + ipji + iz 



\y) 



(6) 



Since the eigenvalues (|lj) are complex, it is important that 
m > max ImAfc in (^-^f. For small p the imaginary parts 
are of order p 2 (second order perturbation theory) so it is 
enough to have m > p 2 in units where the infrared scale 
is of order 1. For large p, m should be made large and 
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then reduced after integration. This will be understood 
throughout. 

Setting A1.2 = A ± A/2 and neglecting the effects of A 
in the averaging in (||) 0], we find that in the flavor sym- 
metric limit, the correlation function in_f^|) relates to the 
'baryonic' pion correlation function |[10fl after a proper 
analytical continuation of the current quark mass 0. 
Specifically, 
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(7) 



where 



l ab C„ B (x,y;z) = 

Tr (S(x, y; z, fi)ij 5 r a S(y, x; z, -fi)i"/r oT b ) ) (8) 



with z — m — iX/2 and E = J dA. For conventional pions 
both propagators in (^) carry ^ with the same sign. 

3. The effects of /i in ([[]) is that of a complex and 
constant 4-vector potential A4 = ifi. It breaks particle- 
antiparticle symmetry much like a vector potential breaks 
particle-particle (antiparticle-antiparticle) symmetry, ft 
acts like a complex electric Aharonov-Bohm effect in the 
particle-antiparticle channel. The particle-antiparticlc 
system breaks apart for fi typically of the order of the 
binding energy (about the pion mass). This phenomenon 
is reminiscent of the destruction of heavy-mesons by 
chromo-electric fields | flT| , and charge or spin density 
waves by transverse electric fields ]l2| , although not iden- 
tical since in our case the 'electric field' is zero. 

In the quenched approximation, the only dependence 
on /i in (||) is that shown in the external propagators. In 
the semi-classical approximation and for three colors (for 
two colors see below) it acts as a 'complex' flux on the 
'diffusons' (particle-antiparticle) 0. For z = m, the long 
paths contributions to (@) are given by 
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(9) 



with Q a = n a 2n/L and Q a = Qa + 2i\ib/± a in V = L 4 . 
The factor 2 in front of A4 reflects on the fact that the 
fluxes add in the 'diffuson'. Using the Gell-Mann Oakes 
Renner (GOR) relation F 2 m 2 r = mS, and the analytical 
continuation m — > m — iX/2, we find 



C nB (x, y; z) w 



JQ-(x-y) _ 



2E 



-iX + 2m + DQ 2 



(10) 



with the diffusion constant D = 2F 2 /Y, 0. Inserting 
( |io| ) into (Q), we observe that the 'diffusion' pole in the 
lower part of the complex plane depends critically on the 
value of the chemical potential /j,. 

In the zero mode approximation n a = or for large 
times t > T crg ~ L 2 /D, the quark return probability is 



p(t,fi) w 6{m^ - 2/i) e -^(»4-4 M 2 )|t| 



(11) 



where we have used E/A = N and g — l/AV, with 
£ = irg, according to the Banks-Casher relation. Here A 
is the mean interlevel spacing between the eigenvalues for 
fi = 0. The occurrence of the step-function theta in ( [Tl] ) 
reflects on the fact that the diffusion pole moves from 
the lower-half to the upper-half of the complex A-plane. 
For t > T crg the quark return probability vanishes for 
H = m^/2 in the quenched approximation. Physically, 
this means that the complex electric flux splits the quark- 
antiquark pair in the quenched approximation a situation 
reminiscent of the magnetic fluxes in type-I superconduc- 
tors |l3) . This result is consistent with current quenched 
lattice simulations H and the results of schematic chiral 
random matrix models for finite [i |l4|-|l6|] . 

We note that in the double scaling limit Dm 2 tH ~ 
mV < 1 and Dfi 2 t H ~ p?V < 1, (|ll|) is about 1 and 
universal. This regime is amenable to a random matrix 
model analysis and signals the onset of a new universality 
for the complex eigenvalues of ([!]). It can be modeled 
using a reduction to 0-dimension. 

4. For two-color QCD the situation is special since in 
this case the Dirac operator possesses an additional sym- 
metry JlTj due to the pseudo-real nature of the SU(2) 
representations. In the diffusive picture of the QCD vac- 
uum this means that the quark return probability re- 
ceives contributions from both 'diffusons' (diagonal) and 
'cooperons' (interference) paths in the semi-classical ap- 
proximation |Q (and references therein). The 'cooper- 
ons' are just the weak-localization contribution to the 
quark return probability ||. In the standard description 
of diffusion they follow from the interference between the 
classical loops traveled in opposite directions (coherent 
backscattering) and reflect on the time-reversal invari- 
ance of the underlying microscopic Hamiltonian. While 
the 'diffusons' sense 2 flux lines, the 'cooperons' are flux- 
blind. A rerun of the above arguments now give 
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(12) 



instead of (||). The first term is the 'diffuson' contribu- 
tion, and the second term the 'cooperon' contribution. 
Inserting ( ^2|) into (Q) yields 

p(t, /i) w e - D <W (9(m n - 2 M ) e +4D ^ + l) (13) 

in the zero mode approximation or t > T cr g = L 2 /D. 
For fi > m^/2 the 'diffuson' contribution (first term) 
drops and we are only left with the 'cooperon' contribu- 
tion which is of order e ~' lra \ t \ . The latter is of order 1 
and universal for /i 2 < m <C 1/V. The transition to the 
'cooperon' phase is simply a transition to the supercon- 
ducting phase in this case. 
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In the universal regime and for small current quark 
masses, the chemical potential is p <C l/VV and small. 
Hence, the complex eigenvalues A& carry an imaginary 
part of order 1/V which is of the order of the microscopic 
level spacing for the 'unperturbed' real parts. If we focus 
on the level-correlations between only the real parts of 
A's in the microscopic limit x = VX ~ 1 we expect a 
mutation in the level correlations from the orthogonal to 
unitary ensemble. The mutation follows a migration of 
part of the quark levels from the real axis to the complex 
plane under the influence of the tiny chemical potential. 
The spectral rigidity £2 (N, p) for the real parts of A's 
can be estimated using semi-classical arguments j?J . The 
result is fl 



E 2 (iV,//)« 9(m„ - 2p) ^ln 



N' 2 



2tt 2 



In 1 



(14) 



where N = E/A > 1, a = D{m\ - Ap 2 )/2A and a = 
2m/ A. The level spacing A = 1/gV is taken to be that 
of the p = limit. The behavior ( [l4| ) can be addressed 
using current quenched lattice Monte-Carlo simulations 
in QCD |,|. 

5. In unquenched QCD, the vacuum supports quark- 
antiquark pairs. The 'baryonic' pion correlations are 
screened by pair creation, rendering the quark-antiquark 
system blind to the complex electric Aharonov-Bohm flux 
(constant A4). As a result, the correlations in (||) are 
primarily that of a quark-antiquark in a vacuum for zero 
nucleon density with p < m^/3 where mjv is the nucleon 
mass (ignoring binding energies). 

At finite nucleon density, the quark return probability 
follows from a pertinent analytical continuation of the 
pion propagator in matter. In a mean-field approxima- 
tion we have plpol 
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(15) 



with Q a = n a 2ir/L in V = L 4 . Here a = (N\qq\N)/T, 
measures the strength of the pion-nucleon sigma term 
relative to the scalar condensate, with 1/a ~ 3po- The 
parameters /3 ~ a and 7 ~ 2a relate to the S-wave pion- 
nucleon scattering lengths [§|,|l9| . The leading density ap- 
proximation follows from the mean-field analysis by keep- 
ing only the leading term in the nucleon density p [0,[l9[ . 
In the space-like regime under consideration there is no 
imaginary contribution to ( |15[) . 

Using the GOR relation F^m 2 = m£, and the ana- 
lytical continuation m — > m — iX/2, we may rewrite ( p~5| ) 
as 



C nB (x,y;z) w 



JQ-(x-y) 



2£(l-ap) 



-iX + 2m + D A Ql + D S Q 2 



(16) 



with temporal and spatial diffusion coefficients 



Da = D 
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1 — ap 



(17) 



Hence D4 ~ D and p independent, while D$ vanishes 
for p ~ I/7 ~ l/(2a) ~ l-5po i n the mean-field approxi- 
mation [pTl. At this point the quark density of states at 
zero virtuality is about jAp) ~ (1 — ap)£ ~ g/2 by the 
Banks-Casher relation ]22[ ]. Using the Kubo- formula we 
conclude that the conductivity vanishes along the spatial 
directions as = DgQ = 0. This is not a metal-insulator 
transition as the conductivity 0-4 = D 4 g ^ is still non- 
zero. It can be regarded as an 'asymmetric' percolation 
transition from d=4 to d=l, with a diffusive quark return 
probability M of the form 



p{t,p) 



-2m\t\ 



(18) 



where we have used E/A* = N ^> 1 and a density 
dependent level spacing A* /A ~ 1/(1 — ap). Here 
E4 = D4/L 2 is the 'temporal' Thouless energy as op- 
posed to Eg = Dg/L 2 the 'spatial' Thouless energy. A 
similar phenomenon takes place at finite temperature [0 . 

At this stage there are two courses of action: The 
conductivity (74 vanishes continuously with the depletion 
of the number of quark states at zero virtuality corre- 
sponding to a vanishing of the quark density of states at 
p ~ 3po- (In fact this is what happens if only the lead- 
ing density approximation were used.) This transition is 
likely of second-order or higher and would be in overall 
agreement with some constituent quark model results [g| . 
Alternatively, it may terminate abruptly for a density 
1.5/5o < P < 3/?o through a d=l percolation transition. 
This is intuitively more appealing if we were to proceed 
from the high-density region backward, and support the 
idea that a nucleon Fermi-surface in 3-dimensions cor- 
respond to an array of 'rods' in d=4 Euclidean space 
(nucleon-worldlines) forcing the conductivity to be es- 
sentially 1-dimensional by Pauli-blocking. Ideas in favor 
of a percolation transition at finite density have been also 
stressed recently by Satz [^3| using different arguments. 

6. We have shown that the disordered properties of 
the QCD ground state are quantitatively altered by a 
finite chemical potential p. In quenched QCD the ef- 
fects of A4 = ip are analogous to that of a complex elec- 
tric Aharonov-Bohm effect, causing the 'baryonic' quark- 
antiquark pair to accumulate 2 flux lines and rupture at 
p = m^/2. This result is in agreement with quenched 
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lattice simulations [|||]. In two-color QCD, the quark- 
quark and antiquark-antiquark pairs are flux-blind. As 
a result, the weak- localization contribution to the quark 
return probability remains unaffected. In the universal 
limit /i 2 < m <C l/V, the quark spectrum in two-color 
QCD exhibits a change in the spectral statistics from an 
orthogonal to unitary ensemble again at p = 

In the unquenched approximation the electric flux is 
screened and the diffusion becomes asymmetric. The 
light quarks take longer time to diffuse along the spa- 
tial directions owing to the presence of a Fermi surface. 
This asymmetry is commensurate with the softening in 
the pion dispersion relation at p ~ po due to pion-nucleon 
S-wave rescattering. As a result, the bulk Ohmic conduc- 
tivity of the system becomes quasi 1-dimensional, with a 
potential for a d=l percolation transition in the density 
range 1.5po < P < 3po- 

Most of our results can be numerically checked by an- 
alyzing the quark return probability in the ergodic and 
diffusive regime at finite chemical potential in lattice 
QCD or in continuum models such as the instanton liquid 
model |@ . 
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